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Получено операторное уравнение для плотности распределения решений системы линейных дифференциальных уравнений с 
полумарковскими коэффициентами, на базе которого выведены зависимости для моментов решений, позволяющие исследо- 
вать устойчивость решения рассматриваемой системы. 

Исследованию устойчивости решений диффе- ной литературе рассматриваются системы диффе- 
ренциальных уравнений со случайными коэффи- ренциальных уравнений, коэффициенты которых 
циентами посвящено много работ [1-5]. В извест- зависят от марковских цепей или марковских не- 
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прерывных процессов и получены условия устой- 
чивости решений как в терминах моментных ура- 
внений, так и в терминах функций Ляпунова. 

В предлагаемой работе рассматривается систе- 
ма линейных дифференциальных уравнений 

^1=А(1,д(1))Х(1), (1) 

аі 

где д(1) - полумарковский конечнозначный про- 
цесс, принимающий конечное число состояний 
в случайные моменты времени 1 0 = О,//,..., 
(?„</</<...). После попадания в состояние Ѳ \ слу- 
чайный процесс д(1) попадает в следующее состоя- 
ние в соответствии с матрицей условных вероятно- 
стей перехода П и независимо от предыдущего со- 
стояния Ѳ к . Описание полумарковских процессов и 
интервально-непрерывных вероятностей приво- 
дится в работах 11,5]. 

Каждому ненулевому элементу л л матрицы 
условных вероятностей перехода 


ставится в соответствие случайная величина Т 1к - вре- 
мя пребывания в состоянии Ѳ к до перехода в состоя- 
ние Ѳ г При этом заданы функции распределения 

РА { ) = Р { Т ,А 1 }’ (к,и = 

Величина Т !к может быть распределена непре- 
рывно или дискретно. Полагаем, что случайная ве- 
личина Т !к неотрицательна и непрерывно распреде- 
лена с плотностью распределения 

/А0 = — (к,з=1,...,п). 

си 

Полагаем также, что полумарковский случай- 
ный процесс д(і) определяется интенсивностями 
перехода из состояния Ѳ к в состояние Ѳ, [1,5] 

чА*) = п ,АА г )’ = 1, ...,«)• (2) 

При этом выполняются условия 

Чл (0 > 0, (? > 0), | ^ (І)сіі = к л ,(з,к = 1,...,гі). 

0 

Пусть д(1) имеет скачки в моменты времени 
И)Ді/ѵч (к~ 0<?і<?2<...). 

Система ур. ( 1 ) распадается на п систем диффе- 
ренциальных уравнений, соответствующих различ- 
ным реализациям случайного процесса д(1) 

^Ш = А к (і)Х к (і), (к = 1,...,п;1>0). (3) 

аі 

Для вывода общих формул предполагаем, что 
для систем (3) известны фундаментальные матрицы 
решений Ы к (і), определяющие решение систем (3) 

X, (0 = Ы к (1)Х к (0), Н к (0) =Е, (к = 1,..., п). 

Если д(1-0)=Ѳ к и д(І к +0)=Ѳ„ то при 1<«1 ]+] си- 
стема ур. (1) принимает вид 

*Ш = А і (і-і ] )Х(і). ( 4 ) 


Будем также предполагать, что в момент / скачка 
процесса д(і) решение системы уравнений (1) имеет 
скачок, определяемый векторным уравнением 


Х(/ ; + 0) = С/Х/ - 0), беі С А * 0, (5, к = 1,..., /7). (5) 


Пусть случайный процесс Х(і),д(1) имеет плот- 
ность распределения 

Ді,Х,д) = ±/ к ((,Х)5(д-Ѳ к ), 

к = 1 


где 5(д) - дельта-функция Дирака. 

Выведем систему уравнений для частных плот- 
ностей / к (1,Х), (к=1,...,п) решения системы линей- 
ных дифференциальных уравнений с полумарков- 
скими коэффициентами (1). Используем вектор 
частных плотностей вероятностей 


Р{і,Х) 


т,х) 


/А А) 


и рассмотрим последовательность векторов Щ,Х) 
(/=0,1,2,...), где / - моменты скачков полумарков- 
ского процесса д(1). В моменты скачков / 
(/=0,1,2,...) вся предыстория случайного процесса 
Х((),д(() "забывается", т.е. не влияет на поведение 
решений системы (1) при />/. Поэтому существует 
стохастический оператор Ь(1)е5$ п1 такой, что 

ГЦ] + 0 = ШАНА (] = 0,1,2,...; I > 0). (6) 

Поскольку все моменты скачков / (/'=0,1,2,...) 
равновероятны, то для простоты изложения, в ка- 
честве начального момента времени возьмем 1 0 =0. 

Система уравнений (6) принимает вид 

Р(1) = Ц1)Г(0), (1> 0) (7) 


или в скалярной форме 

/А,х) = ^ь ь т(0,х), 

5=1 

Ь ь (0 е ,{к = \,...,п-,1> 0). 


Пусть случайный процесс д(1) при /„=0 попадает 
в состояние Ѳ к . При этом выполнены следующие 
условия 

/ДО, X) = 0 (/**), / (0, X) > 0, I / (0, Х)йХ = 1 . 

Е т 


С вероятностью у/ кк (і) процесс д(Г) остается в со- 
стоянии Ѳ к в течение времени 1 > 0 и с вероятностями 
^Дт)т/г в течение временного промежутка [т,г+с!т] 
переходит в состояние Ѳ ; (5=1,... ,я). Кроме этого в 
момент скачка происходит скачок фазового вектора 
Х(Г + 0) = С^Х(Т-О), (5 =1,..., 77). 

Для частных плотностей получим выражение 
Л (Б л = V* (0/ (0, Щ\і)Х) ёеі Щ\І) + 
г п 

+ \ЦчА т ) 1 А ( - 

0 -5=1 

х(0, Х/ДОС/Х) ёеі Х/Дт) Іёеі сАсіх, 
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./; (о у) = } X ч* (44 (* ~ *)Л № 4 1 (ОС^ЛО х (8) 

О 5=1 

хёеІІѴ/(г)|(1е1: С~' к \сІт, (/ Ф к, 1 =1 ,...,«). 

Упростим запись системы ур. (8) с помощью 
введения специальных обозначений. 

Введем в рассмотрение стохастические опера- 
торы К кк е8 + и (к=1,...,п), 4еУ/ (,у,*=1,...,п), опре- 
деляемые формулами 

4Л4 - Л^'(/)У)СІСІ уѴ;'(/), (А: =1,..., п), 

Л/Л) - /(С/У) |сіе1 С/ , 4 к = 1, и). 

Операторы К кк е8\ ь (к=\,...,п) определяют изме- 
нение плотности распределения случайной вектор- 
ной величины Х(Г) при линейном преобразовании 

Х(і) = М к (Т)Х(0), (к = \,...,гі). 

Стохастические операторы К !к е8 + и , (з,к=\,...,п) 
определяют изменение плотности распределения 
при линейном преобразовании 

У = С Л Х, ($,к=\,...,п). 

Систему уравнений (8) можно переписать в виде 

4Л(0Л)=ѵоа(04(О/(0,у) + 

/[Хп 0040-044 (О/ (О ,х)(іт, 

О *=1 

/а./ № У) = |Х 9л 004 С - т)Л Л (т)Х (0, Х)йт, 

0 5=1 

(/ ^ к\к,1 = 1,..., 7г) 
или в операторной форме 

ЛЛ(0,У) = 5 Л Л(ОЛ(0,У) + 

} X у* 004 о -04*4 оол (°. 44 ( 9 ) 

0 »=• 

Введем матрицы, элементами которых являют- 


ся операторы 








Аі(0 4(0 

•• 4,(0 


40 = 


4(0 4(0 

•• 4,(0 
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4(0 О/О 

•• 4(0 





4(0 о 

... 0 



40 

= 

о 4(0 

... 0 

9 




0 0 

-4(0. 



^ІіСО^ІІ #12 (0*5*12 

- 0і,(О4 

40 = 

#21 (0^21 #22 (0*^22 

”• #2 п (0^2 


Ч, 

і(04і ч„і(*)8„2 

- ч„Л 04 


Система ур. (9) выполняется, если выполняют- 
ся операторные уравнения 


Л (0 = <Ѵ/Оа4 (0 + |Хда 004 0-т)4 4 (т)дт, 

0 5=1 

{1,к = \,. ,.,п). 


которые можно записать в матричной форме 

7(0 = ЧД040 + ^ Ь(1 -т)8(т)К(т)сіт , 

0 


где обозначено 


по 


Ѵ4(0 о ... о 
о і// 22 (0 - о 


о 0 ... ц/ пн (1) 


( 10 ) 


Используя равенство (7), можно записать систе- 
му уравнений для вектора частных плотностей Р(Г,Х) 

Р((,Х) = Ч > (іЩі)Р(0,Х) + 

+| Щ - т)8(т)К(т)Р{ 0, Х)с1т. 

0 


Полученный результат сформулируем в виде 
теоремы. 

Теорема 1. Пустъ коэффициенты системы линей- 
ных дифференциальных уравнений (1) зависят от по- 
лумарковского конечнозначного случайного процесса 
ф), который определен заданными интенсивностя- 
ми цф), (з,к=1,...,п). Пустъ между двумя последова- 
тельными скачками случайного процесса ф) при 
!<1<1 к +\ при д(1)=Ѳ , система уравнений (1) совпадает 
с системой (4). Пустъ решения системы уравнений 
(1) имеют скачки вида (5), происходящие одновре- 
менно со скачками процесса д(І). Тогда частные 
плотности /ф,Х), (к=1,...,п) случайного процесса 
(Х(1),д(1)) определяются уравнением 
Р((,Х) = ЦІ)Р( 0, У), 

где оператор Ь(і) удовлетворяет операторному ура- 
внению (10). 

Операторное уравнение (10) можно решать ме- 
тодом последовательных приближений, который 
лишь в исключительных случаях может дать реше- 
ние в замкнутой форме. Преобразуем ур. (10) к бо- 
лее удобному для вывода моментных уравнений 
виду. 

Теорема 2. Решение ур. (10) можно представитъ в 
виде 

Щ) = ЧД 040 + 1 ЧфгЩгМг -т)сіт, ( 11 ) 
0 

где оператор (7(1) удовлетворяет интегральному опе- 
раторному уравнению 

(7(1) = 4040 + 1 8(і-т)Щ-т)іі(т)сіт. ( 12 ) 
0 


Доказательство. Подставим выражение (11) в 
ур. (10). Ур. (10) будет выполнено, если будет спра- 
ведливо равенство 
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І 4 (т)К(т)17(1 -т)</т = 
о 

= Г 4(1 -т)К(і -т)8(т)К(т)с1т + 

о 

( ,-т ^ 

е| I «['(^адС/^-т-# 8(т)К(т)с1т. (13) 


о V о 


Изменим порядок интегрирования в двойном 
интеграле и получим равенства 


( I- X 




8(т)К(т)сіт = 

\ 


| | 4($)К($)Ѵ(і ~Т -$)</$ 
о V о 

| Ч / (5)Л(5)й^5 Г Щ(-Т -5)8(г)К(т)<ІТ 

о V о 

і / 1 

= 1 4(1 -т)К(1 -т)с!т ^11(1 -$)8($)К($)сІ$ 

о Ѵо 

с помощью которых ур. (13) можно записать следу- 
ющим образом 

1 4 (і - т)К(г - т)ІІ (т)сіт = 

О 

= } 4(1 - т)К(і - т)8(т)К(т)с!т + 

О 

+| ѵ Р(г - т)К(і -т)сіт { |[/(г —х)3 ($)К($)сі$ 


Ѵо 


Очевидно, что это уравнение справедливо, если 
Щі) = 8(і)К(і) + } 1/(1 -т)8(т)К(т)<іт. ( 14) 

О 

Будем искать решение операторного уравнения 
(14) в виде 

1/(1) = 8(1)К(1 ) +'\8(1-т)Я(1-т)Ѵ(т)(іт. (15) 

О 

Подставляя ( 1 5) в ур. ( 14), получим уравнение 
1 8(1 -т)К(і -т)Ѵ(т)<іт = 

О 

= '^ 8(1 -т)П(1 ~т)8(т)К(т)с1т + 

О 

1 Г 1-т \ 

+| | 8($)К($)Ѵ(1 - т - $)</$ І8 (т)К(т)сіт . (16) 

О V 0 ) 

Изменяя порядок интегрирования в двойном 
интеграле, получим равенства 

' ( х ~ х Л 

I | 8(х)К($)Ѵ(1-т-з)с1з \8(т)К(т)сіт = 

о V о ) 

' Л-* 

= 1 5(а)К(з)Ж! \Ѵ(1 - т -$)8(тЩт)с1т = 
о ѵ о ) 

= '\8(1-т)К(1-т)сІт[ | Ѵ(і-$)8(з)Я($)сІ$ 


V о 


Используя эти равенства в ур. (16), приходим к 
операторному уравнению для оператора Ѵ(1) 

Ѵ(І) = 8(1)К(1) + 1 Ѵ(1 -т)8(тЩт)с/т. (17) 

О 

Сопоставляя ур. (14) и (17), видим, что можно 
положить П (і)=Ѵ(і). При этом замену (15) можно 
рассматривать как операторное уравнение (11) 

V (!) = 8(1)К(1) + ^(І - т)К(( - т)1/ (т)с/т, 

о 

что и доказывает справедливость теоремы. 

Замечание. Операторное уравнение (14) можно 
записать в виде 

П(і) = 8(1)К(1) + 1 Ѵ(т)8(1 - т)К(1 - т)с!т. 

О 

Сравнивая это уравнение с ур. ( 12), видим, что в 
ур. (12) можно переставить операторы ІІ( т) и 
8(1— т) Щі—т) в подынтегральном выражении. 

Для вывода моментных уравнений умножим 
операторные уравнение (1 1) и (12) справа на вектор 
Р(0,Х) и получим систему уравнений 
Р(1,Х) = 4(і)К(і)Р(0,Х) + 

+| 4(1 -т)К(1 -т)Н(т,Х)сіт, 

Н(1,Х)=8(1)К(1)Р( 0,Х) + 

+ [8(1-т)К(1-т)Н(т,Х)сіт, 

(18) 


где Р(1, X) = Ц1)Р( О, X), Н(1, X) = Ѵ(і)Р( О, X). 
Используя обозначения 


Р(і,Х) = 


'№,х) 


ІЛих) 


, Н(1,Х) = 


Чих) 


К (их) 


можно записать систему ур. (18) в скалярной форме 

мі,х)=ѵ № т к т(о,х)+ 

+\Ѵа ( { ~ Т ) К и ( ( ~ Т Ж ( т . Х )ёт , 

О 

і Ч (і,х)=^ Чь тл,т( о,х)+ 

5=1 

(*,Х)сіт, (20) 


О ■у= 1 


(к = \,...,п). 


Систему ур. (20) для частных плотностей распре- 
деления / к (!,Х), (к=\,...,гі) можно непосредственно 
использовать для вывода моментных уравнений при 
условии, что К кк <=8$, (к=\,...,п), %^8$, (у,*=1,...,и). 

Вектор моментов первого порядка 

М( 1 ) = {Х( 1 ))= \х/( 1 ,Х)сіХ 

Е т 


17 



Известия Томского политехнического университета. 2005. Т. 308. № 4 


можно выразить через частные моменты М к (і), 
(к=1,...,п) первого порядка 

М(Г) = X М к (0, М к {() = | Х/ к (1,Х)с1Х, 0 к = 1,..., л) 

*=! Е т 

так как 

/(і,Х) = ^/ к (і,Х). 

к = 1 

Аналогично матрицу моментов второго порядка 
0(1) = (х(і)Х'(і)) = | XX* / (1,Х)с1Х 

Е т 

можно выразить через матрицы частных моментов 
второго порядка 

т = і ВД, о к (0 = \ хху к (1,Х)СІХ. 

к = ] Е т 

Введем вспомогательные векторы 

Ѵ к (()= | ХИ к (і,Х)сІХ, (к = \,...,гі) 

Е т 

и вспомогательные матрицы 

ж к (о=| хх\ (и х)сіх, (к = і, л). 

Е т 

Умножим систему уравнений (20) на вектор Хи 
проинтегрируем полученные равенства по всему 
пространству Е т . Используем следующие равенства 

| Х/ к (і,Х)сІХ = М к (1),(к =\,...,п), 

Ё т 

| ХК № (і)/ к (0,Х)<ІХ = 

Е т 

= _[ Х/ к (0,Ы- к \і)Х)йеіН- к \і) ах = 

Е т 

= | Х к Ш к (0,У)<іУ = Х к (1)М к (0)ук=\,...,п). 

Е т 

При выводе использовалась замена У=И к (і)Х. 
Аналогично можно получить равенства 

\хК кк (1-т)И к (т,Х)дХ = 

Е т 

= | Хк к (т ,Ы к ' (I -т)Х)&&іЩ 1 (1 -т)сІХ = 

Е ш 

= \ Ы к (1-т)У1ф,Г)с1Г = Х к (1-т)Ѵ к (т), 

Ё т 

(к = \,...,гі), 

/хзд/до ,Х)0Х = 

Е т 

= | х/до,х;'(ОС ь 1 х)ёеі2ѵ; 1 (0|сіеі^ух = 

Е т 

= | ад т: (о, у)йу = ад (о м г (0), 

Ё т 

(к, 5 = 1,..., л), 


\х8 ь Куі-т)к х (т,Х)сіХ = 

Ё, п 

I Хк ч (т,Ы-'( 1 -т)С- к ІХ)йсху'( 1 ~т) ОсІ С ь ' (іх = 

Е т 

= | <ад а - ад (г, уу у = ад (г - ад (г), 

Ё т 

(к,$ =1,...,л). 

Окончательно приходим к системе интеграль- 
ных уравнений 

м к (і)=ѵ № т к тл о)+ 

+}ѵ / й( ; -' г ) ЛГ г( ; -' г )^( т ) Л 

о 

Ѵ к (і) = Хсі ь (і)С ь ху(Щ(0) + 

5=1 

} х ^ г ) с ьЛ О-г)0(гу/г, (21) 

О »=1 

(А: = 1 ,..., й), 

которые определяют частные моменты первого по- 
рядка /Ц,(/), (А=1,...,л). 

Аналогично находим систему матричных инте- 
гральных уравнений для матриц частных моментов 
второго порядка Д .(/), (к=1,...,п) 

о к (()=ѵ № № к ((Щ(т к Хі)+ 

+{ Ѵкк (1-тЩ(1-тЩ(т)Х к (1-т)4т, 

0 

ъ (о = Ё 9ь (ос„ к т (о ж ик + 

5 = 1 

+} Ё ( 1 ~ Т ) С ь Д (І-'С)Щ (т) N‘(1 - г)Д* с/г, (22) 
0 ®=1 

(Лг = 1,..., л). 

Полученный результат сформулируем в виде 
теоремы. 

Теорема 3. Пусть выполнены условия теоремы 1. 
Тогда математическое ожидание 

М({) = {Х(())=Ё М Л*) 

к=\ 

случайного решения системы ( 1 ) определяется систе- 
мой интегральных уравнений (21), а матрица вторых 
начальных моментов 

о( *) = {х«)хХо) = Ё в Л) 

к = 1 

определяется системой интегральных уравнений (22). 

Замечание. Если в системе ур. (1) отсутствуют 
скачки решений (5), то в формулах (21) и (22) сле- 
дует положить С к =Е, (к,х=\,...,п), где Е - единич- 
ная матрица. 
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